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materials,  that  there  exists  a remarkable  universality  of 
dielectric  response  behavior  regardless  of  physical  structure, 
types  of  bonding,  chemical  type,  polarizing  species  and  geo- 
metrical configurations.  This  strongly  suggests  that  there 
should  exist  a correspondingly  universal  mechanism  of  dielectric 
polarization  in  condensed  matter.  The  present  work  proposes 
such  a universal  mechanism  associated  with  the  existence  of 
, some  ubiquitous  very  low  energy  excitations  in  the  system.  These 
j excitations  exhibit  an  infrared  divergent-like  response  to 
transitions  of  the  polarizing  species  induced  by  a time-varying 
electric  field  in  the  dielectric  and  give  rise  to  the  universal 
dielectric  response. 
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FREQUENCY  DEPENDENCE  OF  DIELECTRIC  LOSS 
IN  CONDENSED  MATTER 


I . INTRODUCTION 

The  dielectric  response  of  solids  and  liquids  has  been  the  subject 
of  intense  investigation  over  a long  period  of  time  extending  to  this 
date,  and  pursued  by  physicists,  chemists  and  engineers  alike.  A 
detailed  survey  of  the  dielectric  properties  of  a wide  range  of  solids 
has  been  given  recently  by  Jonscher.^  It  was  observed^  that  the  die- 
lectric response  functions  in  frequency  or  in  time  depart  strongly  from 
the  Debye  response  for  a large  number  of  essentially  dissimilar  materials 
and  fall  into  a remarkably  common  or  "universal"  pattern.  In  particular, 
the  frequency  dependence  of  dielectric  loss  follows  the  empirical  law 

X"  (a))  « U)n  1 with  0 < n < 1 (1) 

extending  over  several  decades  of  frequency  from  low  audio  and  sub-audio 
to  u)/2n  '/'lO  Hz.  For  some  dielectrics,  a broad  loss  peak  may  be  found 
at  lower  frequencies.  Genuine  Debye  behavior  with  the  complex  suscepti- 
bility given  by  x(w)  “ (1  +i(i)T)  * is  seldom  observed  in  solids. 

Examples  of  the  materials  that  obey  the  empirical  law  (Eq.  (1))  include 
inorganic  ceramics;  ionic  conductors;  polymeric  materials,  inorganic 
crystalline  and  amorphous  materials  including  glasses,  insulating  or 
semiconducting;  and  organic  and  biological  systems.  By  way  of  these 
examples  we  see  that  the  frequency  response  (1)  is  similar  for  systems 
with  permanent  dipoles  and  with  hopping  charge  carriers  of  electronic  or 
ionic  nature.  It  is  valid  in  covalent,  ionic  and  molecular  solids,  in 
single  crystals,  polycrystalline  and  amorphous  structures;  hence  the 

behavior  (1)  is  apparently  independent  of  the  particulars  of  the 

• . . 9 

material.  At  higher  frequencies,  10  Hz  and  up,  quantum  effects  involv- 
ing lattice  mode  excitations  and/or  electronic  excitations  become 
prominent  and,  as  is  well  known,  the  response  then  differs  from 
material  to  material,  and  as  such  will  not  be  of  interest  to  us  in  the 
Note:  Manuscript  submitted  August  16,  1978. 
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present  context. 

The  various  types  of  dielectric  response  are  summarized  in  Fig.  1. 
We  note  the  virtual  absence  of  the  pure  Debye  response;  and  the  valid- 
ity of  the  universal  law  of  dielectric  response,  (1),  in  a remarkably 
wide  range  o physical  and  chemical  situations,  and  over  a very  wide 
range  of  frequencies.  In  some  types  of  dielectrics  the  universal 
response  (1)  is  followed  at  low  frequencies  by  a loss  peak  referred  to 
as  a and  0 peaks,  or  by  another  universal  response  with  n typically 
between  0.1  and  0.3. 

It  is  this  state  of  affairs  that  has  motivated  us  to  seek  a 
renewed  understanding  of  these  phenomena  in  terms  of  a common  or 
"universal"  characteristic  across  the  entire  spectrum  of  materials  and 
to  associate  such  a characteristic  with  some  physically  simple  and 
"elementary"  principles  or  properties.  In  the  next  section  we  shall 
present  several  elementary  principles  which  when  combined  enable  a 
derivation  of  the  "universal"  law  (1)  regardless  of  the  physical, 
chemical  and  geometrical  properties  of  the  solids,  and  also  regardless 
of  the  nature  of  the  electrically  active  species  responsible  for  polari- 
zation, whether  dipoles,  electrons  or  ions.  Then  in  Section  III  we 
discuss  several  examples  of  low  energy  excitations  expected  in  a host  of 
systems  that  satisfy  these  elementary  principles.  In  Section  IV  we 
derive  the  universal  response  (1)  and  the  possible  presence  of  a loss 
peak  at  lower  frequencies.  Finally,  in  Section  V we  make  some  con- 
cluding remarks. 


II.  INFRARED  DIVERGENCE  AND  THE  "UNIVERSAL"  LAW  y"  (ui)  « CJ11"1 

Infrared  divergence  phenomena,  although  not  commonly  observed  in 

physics,  have  been  seen  in  several  instances.  The  most  well-known  case 

is  in  quantum  electrodynamics  where  the  infrared  divergence  manifests 
. . . 2 

itself  in  a Bremsstrahlung  experiment  of  a fast  charged  particle.  In 

3 4 

the  realm  of  solid  state  physics  ’ an  example  of  infrared  divergence  i3 
thought  to  be  provided  by  the  peculiar  shape  of  X-ray  absorption  edges 
of  metals,^  and  the  "orthogonality  catastrophe"  for  an  impurity  poten- 
tial inserted  in  a Fermi  gas.^  These  examples  are  by  no  means  exhaus- 
tive but  the  subjects  they  cover  demonstrate  that  infrared  divergence  is 

3 4 

not  uncommon.  Excellent  reviews  on  the  subject  are  available.  ’ 

The  features  common  to  systems  exhibiting  the  infrared  divergence 
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phenomenon  are  (a)  the  sudden  application  of  a potential,  or  a sudden 
change  of  the  potential  or  the  Hamiltonian;  and  (b)  availability  of  low- 
energy  excitations  of  the  system  and  its  response  to  the  sudden  potential 

change  dominated  by  the  emissions  of  these  low-energy  excitations.  In 

. . . 7-9 

the  time  domain  the  phenomenon  is  the  transient  response  of  the 

system  to  that  abrupt  change  of  potential.  Infrared  divergence  occurs 
whenever  the  suddenly  switched  on  potential  V excites  some  low  energy 
excitations,  with  density  of  states  Nos)  for  excitation  energy  E,  which 
is  such  that  V^(E)  |\J(E)  « E.  In  this  instance  there  is  an  increasingly 
high  probability  of  exciting  decreasingly  small  energy  excitations  and 
this  causes  a power  law  divergence  of  the  response  in  the  frequency 
domain.  In  the  X-ray  edge  problem  in  metals  an  X-ray  photon  when 
absorbed,  suddenly  switches  on  a hole-core  potential  V for  the  conduc- 
tion electrons.  The  low-energy  excitations  here  are  the  electron-hole 
pairs . 

In  the  later  sections  we  shall  argue  that  within  a broad  classifi- 
cation of  dielectrics,  according  to  a scheme  to  be  outlined,  there  exist 
states  which,  for  convenience,  we  shall  refer  to  as  correlated  states. 
These  correlated  states  have  a smooth  and  continuous  density  of  states 
N(E),  which  we  take  as  constant  N.  Low  energy  excitation  of  the  cor- 
related states  with  excitation  energy  E consists  of  removing  an 
"occupied"  state  to  an  "unoccupied"  state  and  is  the  analogue  of  the 
electron-hole  pair  excitation  in  the  X-ray  edge  problem. 

The  charged  particles  or  dipoles  responsible  for  polarization  in 
the  dielectrics  undergo  quantum  transitions,  including  changes  in  their 
positions/orientations,  between  preferred  states  in  an  abrupt  manner  by 
hopping  or  jumping  movements  such  that  the  time  1/v  taken  by  the  actual 
transition  is  negligible  in  comparison  with  both  (i)  the  time  spent  on 
average  in  the  respective  preferred  states,  and  (ii)  the  time  charac- 
teristic of  the  low  energy  excitation  of  the  correlated  states.  The 
condition  (i)  is  invariably  satisfied  in  solid  dielectrics.  That  con- 
dition (ii)  is  also  satisfied  will  become  clearer  after  we  have  con- 
sidered the  nature  of  the  correlated  states. 

Due  to  the  charged  particle  (dipole)  transition  a potential  is 
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suddenly  switched  on  which  acts  on  the  correlated  states.  The  low- 
frequency  response  of  the  dielectric  to  this  potential  involves  the  emis- 
sion of  low-energy  excitations  of  the  correlated  states.  We  shall  argue 
that  the  low  energy  excitations  of  these  correlated  states  have  a density 
of  states  |\|(E)  x E,  and  that  the  potential  change  V has  little  or  no  E 
dependence.  It  follows  that  the  conditions  for  an  infrared  divergent 
dielectric  response  of  the  correlated  states  are  satisfied.  The  mean 
number  n of  correlated  state  excitations  is  then  n « bV  N J c EdE/E  , 
which  diverges  logarithmically,  where  Ec  is  the  upper  "cut-off"  of  the 
correlated  state  excitation  energy.  The  Fourier  transform  to  the  time 
domain  of  the  "universal"  relation  (1)  is  i(t)  « t n,  i.e.,  the  widely 
observed  Curie-von  Schweidler  law*  of  depolarization.  It  is  interesting 
to  note  that  the  infrared  divergence  problem  when  considered  in  the 
time  domain  as  a transient  response  problem  does  lead  to  the  time 
decay  of  the  response  function  for  large  times  as  S(t)  « t n.  The 
derivation  of  the  complete  dielectric  response  will  be  deferred  to 
Section  IV,  after  we  have  discussed  the  correlated  states  in  a broad 
classification  of  dielectrics  in  the  next  section. 

III.  Correlated  States 

In  the  preceding  section  we  connected  the  "universal  law,"  Eq.  (1) 
to  an  infrared  divergent  response  of  correlated  states.  For  this  inter- 
pretation to  follow  it  is  necessary  for  such  states  to  be  prevalent  in 
dielectrics  and  satisfy  two  rather  stringent  but  interrelated  criterion: 
namely  (i)  the  characteristic  response  time  for  the  correlated  states  is 
long  in  comparison  to  the  switching  of  the  Hamiltonian  and  (ii)  these 
low  energy  excitations  are  not  in  general  thermodynamically  accessible 
at  the  moderate  temperatures;  i.e.,  several  100°K,  since  otherwise  we 
expect  significant  temperature  smearing  of  the  effect.  Below  we  discuss 
several  examples  of  such  correlated  states  which  can  reasonably  be 
expected  to  be  present  in  many  dielectrics.  One  should  note  from  the 
onset  that  although  these  examples  are  quite  general,  they  are  by  no 
means  exhaustive. 
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(a)  Dielectrics  with  Electron  Pairing  Interactions  and  States 

Anderson^  has  recently  proposed  a model  of  amorphous  semicon- 
ductors which  emphasizes  the  role  of  electron  pairing  interactions  in 
these  systems.  Anderson's  ideas  should  apply  to  a large  number  of  dia- 
magnetic dielectrics.  The  concept  of  strong  local  electron  pairing 
interactions  arises  from  the  observation  that  by  and  large  if  a par- 
ticular electronic  state  is  singly  occupied,  the  atom  or  atoms  prin- 
cipally associated  with  this  state  will  adjust  their  positions  in  such 
a way  as  to  lower  the  energy  of  this  state  relative  to  its  value  if  it 
were  constrained  to  be  unoccupied.  This  effect  can  produce  a self- 
trapping of  the  electron  in  a manner  related  to  the  formation  of  a 
small  polaron;  however,  it  also  makes  it  favorable  for  the  spin-mate  of 
the  occupied  state  to  be  occupied  and  in  turn  the  two  electrons  can  be 
more  strongly  self-trapped  in  a correlated  way.  The  resultant  quasi- 
particle (the  two  electrons  and  the  concomitant  lattice  distortion)  has 
been  termed^  for  obvious  reasons  a bipolaron.  Of  course  the  formation 
of  such  a quasi-particle  is  resisted  by  the  mutual  coulomb  repulsion  of 
the  constituent  electrons,  but  only  in  comparatively  few  systems,  such 
as  the  transition  metal  oxides  or  impurity  bands  in  crystalline  semi- 
conductors does  this  repulsive  interaction  dominate,  giving  rise  to  e.g. 
spin  density  fluctuations. 

In  order  to  develop  further  the  idea  of  local  distortion  mediated 
effective  electron-electron  attractive  interaction  Anderson  has  employed 
an  effective  negative  U Hubbard-like  term  Un. ,n. , to  model  the  effect 
where  tL  is  the  number  operator  for  an  electron  of  spin  o in  a state 
centered  at  the  "site"  i;  |ia>.  One  should  keep  in  mind  that  i could 
well  index  the  up  and  down  spin  states  associated  with  a group  of  atoms 
and  not  just  a single  one.  Let  us  describe  a group  of  such  centers  in 
contact  with  one  another  as  well  as  alternate  states  of  the  system 
(those  with  = 0)  by  the  simplified  Hamiltonian: 

(2) 


+ .?_R. .a*_a ._  + ? U.n..n.,  , 

l i it  i+  ’ 


H = 2 E.n.  _ • . ....  .-. 

1.0  l lG  1JG  lj  lO  jo 


where  a.  , a.  create  and  annihilate  electrons  of  spin  o in  the  state 
i o ia  r 

|io>>  and  we  take  R^  as  R if  i,j  are  nearest  neighbors  and  zero  other- 
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wise . 


The  parameters  {e.}  and  {u^}  are  considered  as  random  variables 
obeying  the  joint  probability  distribution  P(E^,  !h)  which  for  the  time 
being  is  left  unspecified.  This  model  can  be  made  to  mimic  many  dif- 
ferent situations  depending  on  the  choice  of  P(E.,U.). 

^ll1 

To  obtain  results  from  (2)  we  have  developed  a generalized  mean 

field-like  method  which  entails  linearizing  the  many-body  terms  U.n.in., 
XX  XX  XX  llTlV 

as  t^n^n^  ^ 2CTUini-ania_Uini+ni-l-where  U^n^n.^  is  present  to  prevent 
double  counting  of  the  interaction  and  the  conditionally  averaged 
number  of  spin  a electrons  at  the  site  i is  given  by  the  relation 

n.  ■ - — JdEf(E)g.  (E+)  (3) 

with  f the  fermi  function  and  E+  denotes  the  lim  (E+is).  Equation 

S"T0+ 

(3)  provides  a set  of  generalized  Hartree  Fock-like  self-consistent 

. . . s 

relations  determining  the  parameters  n^  since  the  Green's  functions 

g^a  entering  these  formulae  are  defined  as  g^c  = <ij| (z-Hg^)  ^|oi> 

where  H . = ? (E.+U.n.  /T)r..  + .S_R.  .aL,a._  and  hence  depends  on  n.  . 
el  iO  i i 1-0  ia  ijO  ij  ia  jo  r 10 

For  convenience  (3)  can  be  recast  as  n.  = -(Im/u)  J’dEf(E)/(E-E.-U.'fi.  - 

10  1 l 10 

A ^),  where  A^  is  the  usual  self-energy  and  is  a function  of  {n^} , E 
and  R. 

To  solve  these  self-consistent  conditions  we  employ  the  long 

12 

established  coherent  potential  approximateion  (CPA)  to  obtain  the  self 

energies  A..  This  entails  defining  an  effective  medium  characterized 

by  a single  potential  Zq  which  is  energy  dependent  and  can  be  complex, 

in  such  a way  that  the  G(E-Z  ) a <g-^(E)>,  where  the  brackets  < > 

o °iO  av:  av. 

here  and  henceforth  denote  an  average  over  the  random  variables  entering 
and  G is  the  Green's  function  obtained  by  replacing  the  site 
diagonal  random  potentials  of  by  Zq  at  each  site  i.  The  CPA  is 
exact  in  both  strong  and  weak  scattering  (virtual  crystal)  limits  and 
hence  provides  an  interpolation  scheme  for  treating  the  intermediate 
cases.  The  use  of  this  method  greatly  simplifies  the  computations  since 
its  functional  form  of  A^  can  be  easily  found  using  established  tech- 
niques once  the  "lattice"  structure  is  specified.  For  example,  if  we 
assume  a simple  chain  then  A ■ (E-Z^)  -/{(e-Zq)^-4R^( . The  CPA  equation 
defining  Zq  and  hence  A can  be  written  explicitly  for  the  present  model 
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, 


We  now  have  as  inputs  into  the  formalism  some  specified  tempera- 
ture, T,  and  number  of  electrons  in  the  band,  Ng^,  as  well  as  particular 
functional  forms  for  G(E)  and  P(E^,IK).  The  calculation  then  proceeds 

as  follows:  First  we  assume  the  function  n.  (E.,U.)  and  then  solve  the 

ia  11 

CPA  equations  using  a modified  Newton-Raphson  technique  to  obtain  £q(E) 
and  hence  A(Zq).  There  are  presumably  many  £q(E)  satisfying  the  CPA  con- 
dition for  a particular  E;  however,  the  correct  branch  of  E (E)  goes  as 

X ...  . 

<E£+U^n^>  at  large  energies,  and  this  solution  can  be  analytically 

extended  into  the  region  of  interest  (energies  within  the  band)  by  use 
of  the  Newton-Raphson  method.  In  doing  this  we  have  found  it  most  effi- 
cient to  follow  a path  in  the  complex  plane  slightly  above  the  real  axis 
(E  + if);  H * .05)  and  then  take  the  limit  ’T+0+  numerically.  Having  so 
determined  ) we  can  then  find  the  chemical  potential  y of  the 

system  from  the  usual  relation  1^=  (2/ff )Im/dE/{ (e~^ ^^+1)(E+-Eo-A(Eq))} , 
where  in  obtaining  this  condition,  we  have  employed  the  CPA  equation. 

Note  the  CPA  determined  satisfies  the  important  sum  rule 
J*dE/(E-Z  -A(£  ))  = 1.  Having  £ (E+)  and  y,  we  then  calculate  a new 
function  n/^E^jU^)  from  (3)  and  this  procedure  is  iterated  until  self- 

I 

consistency  is  established,  i.e.  n.  „(E.,U.)  = n-  (E-,U.).  If  a con- 
tinuous  probability  distribution  is  assumed  for  the  random  variables, 
it  is  of  course  not  numerically  feasible  to  establish  self-consistency 
at  each  point  in  (E.,  U-)  space.  In  these  instances  we  establish  self- 
consistence  at  a grid  of  points  assuming  that  n^tE^,!!^)  can  be 
adequately  represented  for  intermediate  values  by  trapezoidal  interpola- 
tion. We  have  found  for  simple  continuous  probability  distributions 
that  this  procedure  converges  very  nicely  (well  within  the  realm  of 
numerical  feasibility)  as  the  number  of  points  in  the  grid  is  increased. 
Note  that  usually  more  than  one  self-consistent  solution  exists  and  this 
will  prove  important  to  our  subsequent  development. 

Before  detailing  some  of  the  examples  that  we  have  treated,  it  is 
convenient  to  backtrack  somewhat  and  draw  a relationship  between  elec- 
tron pairing  interactions  and  covalency  and  in  so  doing  motivate  these 
cases  and  further  stress  the  generality  of  the  pairing  ideas. 
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As  a prototype  consider  a simple  dangling  bond  such  as  one  associ- 
ated with  an  Si  atom  which  is  bonded  to  three  neighboring  silicons 
leaving  a dangling  hybrid.  If  we  denote  by  x^  the  displacement  of  this 
atom  from  where  it  would  sit  if  the  dangling  hybrid  were  constrained  to 

be  singly  occupied  with  energy  E,  then  a Hamiltonian  partially 

n 13 

describing  the  energetics  of  this  atom  is  H^=  ^E^n^  -X^x^Cn^+n^l) 

+ c^xj^/2,  where  n^a  is  the  number  operator  for  electrons  of  spin  a in 

the  dangling  hybrid  orbital  |ha>.  The  last  term  entering  is  a back- 

bond  stretching  energy  and  the  second  is  the  so-called  dehybridization 
14 

energy.  If  we  have  a group  of  such  nonbonded  states  interacting  with 

one  another  then  a simplified  Hamiltonian  describing  the  situation  is 

H = ? E^n.  ?X.x. (n.  +n. .-1)  + Ec.x?/2  + .E  R. .a.  a._  We  can  now  view 

ict  i ict  i i i i+  i|  lit  ijCT  ij  i a j a 

the  displacements  x^  as  parameters  entering  the  Hamiltonian  to  be  deter- 
mined self-consistently  by  requiring  the  free  energy  of  the  system  to 
be  stationary  with  respect  to  their  variations.  This  results  in  a set 
of  self-consistent  conditions  which  can  be  used  to  eliminate  the 


parameters  x. 


It  is  then  a simple  matter  to  show  that  the  resultant 


Hamiltonian  is  essentially  similar  to  the  negative  U model  (2)  within 
the  context  of  our  mean  field  approximation  if  we  make  the  identifica- 
tion:  2X^/c.  ■*'U^,E^">,E^+U^/2.  Thus  we  expect  the  negative  U model  to 

incorporate  the  behavior  of  a simple  nonbonded  orbital  since  the 
analysis  of  course  is  not  limited  to  only  the  Si  dangling  hybrid  but 
applies  whenever  one  has  a dangling  bond  associated  with  covalent  back- 
bonds  and  hence  is  quite  general.  Furthermore,  in  amorphous  partially 
covalent  materials  because  of  the  more  localized  nature  of  the  wave 
function  in  comparison  to  the  crystalline  case  one  expects  similar 
effects  to  arise  from  states  with  origin  in  the  conducting  band. 

Indeed  in  an  amorphous  material  one  would  expect  after  Anderson^ 
that  the  one  electron  potential  E.  obeys  a continuous  probability 
distribution  P(E^)  spanning  the  forbidden  gap.  To  model  the  resultant 
situation  we  have  solved  the  self-consistent  equation  (3)  assuming 
P(E^,U^)  - 5(IK-Uq)W(E^)  where  W(E^)/B  = 1/2  for  -1  <_  E^/B  < 1 and  zero 
otherwise  and  Uq/B  = -3;  B is  unperturbed  half  bandwidth.  Also  we 
assumed  that  T=0,  Ne^=l  and  employed  as  an  unperturbed  Green's  function, 
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G,  appropriate  for  a Cayley  tree  of  coordination  number  six.  In  Fig.  2 

we  exhibit  the  numerically  determined  lowest  energy  state  of  the 

system  (solid  line),  as  well  as  another  self-consistent  solution  (dashed  line) 

which  represents  a low-lying  excitation  of  the  system.  The  two  solutions 

essentially  differ  from  one  another  by  the  transfer  of  electrons  from 

one  group  of  pairing  centers  to  another,  and  in  this  way,  although  there 

is  a large  gap  in  the  one-electron  spectrum  very  low  lying  excitations 

can  be  achieved  leading  to  a gapless  pair  state  spectrum.  To  understand 

this  behavior  further  consider  two  isolated  pairing  centers  labeled  i,j 

in  competition  with  one  another  for  two  electrons.  Then  if  U.  * U.  it 

1 J 

is  not  the  magnitude  of  U that  determines  the  occupancy  but  rather  E^, 

E j . For  example  if  E.  <Ej  then  the  site  labeled  i is  doubly  occupied 
and  that  labeled  j is  doubly  empty  in  the  ground  state.  Thus,  although 
the  one  electron  states  lie  at  E^  + and  E j , and  are  hence  usually 
well  separated  in  energy  (^.lev),  excitations  of  the  system  that  require 
only  energy  E^-E^  which  becomes  vanishingly  small  as  E.^Ej  can  be 
achieved  by  removing  the  electron  pair  from  the  site  i to  the  site  j. 

In  the  case  R.  .-*-0,  the  density  of  pair  state  excitations  is  then 
J'  P(E^  +11/2)  which  is  continuous  and  slowly  varying  around  E^  and  as 
we  have  seen  a similar  picture  also  applies  if  we  assume  some  coupling 
between  the  pairing  centers. 

Such  a smooth  distribution  of  self-trapped  pair  state  excitations 

is  expected  to  have  a character  sufficient  to  produce  an  infrared 

divergence  at  very  low  temperature.  This  is  so  because  the  density  of 

2 

states  of  low  energy  pair  state  excitations  with  energy  E is  -''N  (E^)E 
where  N(E^)  is  the  density  of  pair  states  at  E^.  Furthermore,  the 
matrix  elements  V. . of  the  potential  change  V abruptly  switched  on  by 
the  polarizing  species  should  on  the  whole  be  independent  of  Ej-E^,  and 
the  fact  that  the  pair  states  are  strongly  self-trapped  implies  that 
their  response  time  can  be  much  longer  than  the  time  characteristic  of 


the  hopping  or  reorientation  of  the  charge  species.  Thus  all  conditions 
for  an  infrared  divergent  dielectric  response  are  apparently  satisfied. 
At  elevated  temperatures  the  infrared  divergence  will  survive  provided 
the  low  lying  pair  state  excitations  are  thermodynamically  inaccessible, 
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but  this  is  a natural  consequence  for  large  enough  (|uJ/kT»l) 
which  should  produce  an  energy  barrier  between  pairing  states  suffi- 
cient to  effectively  presevent  thermally  assisted  tunneling  in  physical 
systems . 

So  far  we  have  considered  only  pair  state  excitations  that  involve 

essentially  transferring  a pair  of  electrons  from  one  pairing  center  to 

another.  However,  there  is  another  (not  completely  orthogonal)  class 

of  low-lying  excitations  associated  with  pairing  interactions  that  arise 

from  breaking  the  pair  and  "placing"  the  electrons  in  states  associated 

with  non-pairing  sites.  For  example,  in  the  case  of  metal-semiconductor 

(Schottky)  contacts  one  could  envision  transferring  the  electrons  from 

pairing  centers  in  the  semiconductor  (say,  nonbonded  orbitals)  to  the 

Fermi  sea.  A particular  example  derived  from  our  general  model  is  shown 

in  Fig.  3 where  we  display  two  different  self-consistent  solutions  to 

(3)  obtained  by  using  the  previously  detailed  formalism.  We  have  taken 

as  inputs  in  calculating  these  examples:  T=*0,  N|?. 3 and  G appropriate  for 

a Cayley  lattice  of  coodination  number  six.  The  form  of  P(E^,U£)  is 

chosen  so  that  P(E.,U.)  * x<5(U.)  6(E.-E  ) + (l-x)<5(U.-U  )W(E.)  where 
11  1 10  101 

W(EiVB  = 5 (B  the  unperturbed  half  bandwidth)  for  .1  <_  E^/B  .3  and 

zero  otherwise  and  x (the  concentration  of  pairing  centers)  = .1  with 

U /B  = -1.6.  The  solid  line  of  Fig.  3 is  within  our  formalism  the 
o 

density  of  states  corresponding  to  the  numerically  determined  lowest 
energy  state  of  the  system  while  the  dashed  line  represents  a low-lying 
self-consistently  obtained  excited  state;  a fact  that  we  have  verified 
directly  by  comparing  the  energies  of  the  two  cases.  These  two  solu- 
tions differ  from  one  another  by  the  transfer  of  electrons  from  the 
pairing  centers  (which  when  occupied  in  this  example  form  a band  of 
states  j'Uq/2  below  E^  as  shown  in  Fig.  3)  to  the  main  band  with  the 
unoccupied  pairing  levels  now  appearing  ^Uq/2  above  E^.  This  is  exem- 
plified in  Fig.  3 in  going  from  the  ground  to  excited  state  by  the 
slight  increase  in  'E^',  as  well  as  the  decrease  in  the  measure  of  the 
pair  band  below  the  main  band  edge  and  a concurrent  increase  in  the 
measure  of  the  main  band. 

Further  insight  into  this  behavior  can  be  gained  by  considering  a 
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single  pairing  impurity  in  a tight  binding  lattice.  The  situation  can 

t 


be  described  by  the  model  Hamiltonian  H 


4 E n.  ♦ .4  R.  . a'  a.  + -. 

ia  o ia  ija  ij  ia  ja  | 

U-n.,n..+  ^E.n.  , which  represents  of  course  a special  case  of  (2). 

J J+  j+  a J ja 

Approximating  IK  n^n^  , as  before  the  free  energy  of  the  system  as  a 
function  of  n.  can  be  expressed  as 

(Un.+E.)(3G.(£+)/3f )d£ 


F 


K - 


u£?  - (2/u)Im  J3aif(E)  rE  — V *4^“ 

J ” l-(Un .+E . )G . (E+) 


(4) 


J J J 


where  K is  a constant  independent  of  n.  and  we  have  used  the  up-down 

J j.  » 

spin  symmetry  present  for  Uj  < 0 to  replace  n^  by  n j . The  free  energy, 
F,  possesses  a double  minimum  as  a function  of  n^  when  Ej  + U./2  lies 
in  the  vicinity  of  E^  and  U./R  is  »1,  as  is  shown  in  Fig.  4.  In 
arriving  at  these  results  we  have  chosen  for  simplicity  a rectangular 
density  of  states  of  half-bandwidth  B to  model  the  main  band  i.e. 

G^(Z)  = -(l/2B)ln  |(Z-B)/^Z+B)|  and  neglected  temperature  effects  which 
are  unimportant  at  moderate  temperatures  for  physically  expected  U.; 

. > . , X ^ 

i.e.  U.  J'  .1  ev.  At  the  minima  n^  satisfies  the  appropriate  form  of 

Eq.  (3)  and  hence  represents  the  self-consistently  obtained  average 

number  of  electrons  of  one  spin  species  at  the  site  j.  The  two  minima 

hence  correspond  to  distinctly  different  occupancy  of  the  pairing  center 
'J'  s 

since  in  one  case  n.  s 0 and  the  other  n.  s 1.  That  is  on  one  hand 
J J 

almost  two  electrons  occupy  the  pairing  levels  which  lie  approximately 
at  (Ej-Uj)  while  on  the  other  the  pairing  center  is  effectively  unoc- 
cupied and  its  associated  states  lie  at  ^ j • One  can  easily  show  that 
the  two  minima  are  separated  for  large  IK  by  ^ J E j + IK/2  - E^|  and  hence 
such  a negative  U center  can  give  rise  to  a low  lying  excitation  of  the 
system  if  its  characteristic  parameters  are  such  that  (Ej+lK/2)  x E^. 
This  is  of  course  consistent  with  the  previously  obtained  results  sum- 
marized in  Fig.  4 and  is  just  a rather  more  specific  case. 

Although  within  the  context  of  the  present  mean-field  like  approxi- 
mation we  cannot  make  a new  linear  combination  of  the  two  states 
represented  by  the  essentially  degenerate  generalized  Hartree-Fock  self- 
consistent.  solutions  (associated  with  the  minima  of  Fig.  4)  that  reduces 
further  the  energy  of  the  system  (there  is  an  orthogonality  theorem^) 
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such  an  effect  of  course  physically  exists.  The  resultant  intrinsic 
matrix  element  connecting  these  states  should  itself  be  a random 
variable  because  of  the  different  allowed  choices  of  E.,  0.  sufficient 
to  produce  the  same  degree  of  degeneracy.  Such  being  the  case,  one 
expects  the  density  of  the  very  low-lying  excitations  at  a particular 
E to  behave  as  E and  contribute  an  infrared  divergent  dielectric 
response  (1).  We  will  postpone  details  of  this  argument  until  the  next 
subsection  where  we  are  confronted  with  an  analogous  problem  in  terms 
of  tunneling  modes. 

The  essentials  of  the  present  low-lying  pair  state  picture  should 
not  be  smeared  out  at  reasonable  temperatures  since  although  e.g.  the 
details  of  Fig.  3 may  be  somewhat  different  at  different  temperatures 
one  still  finds  a double  minimum  in  F(n^)  and  the  corresponding  low-lying 
excitations,  and  once  again  for  large  thermally  assisted  tunneling 
from  one  state  to  the  other  would  not  be  expected  to  occur. 

Thus,  we  have  detailed  two  rather  general  examples  which  illustrate 
how  electron  pairing  interactions  can  provide  correlated  states  with 
characteristics  sufficient  to  produce  the  "universal  law";  Eq.  (1).  Of 
course  these  examples  are  not  completely  unrelated,  and  one  expects  in 
many  systems  that  both  types  of  excitations  are  simultaneously  present 
and  operative  in  producing  an  u)n  * behavior  of  x" • Although  we  have 
phrased  our  discussion  primarily  in  terms  of  amorphous  systems  where  one 
expects  an  appreciable  number  of  weaker/stronger  bonds,  lone  pairs,  etc., 
to  be  present  giving  rise  to  the  gapless  correlated  states,  it  is  also 
reasonable  to  expect  that  such  low-lying  excitations  occur  and  are 
important  in  more  nearly  crystalline  covalent  solids  since  the  remaining 
pairing  centers  in  these  materials  could  partially  pin  the  Fermi  level 
in  their  vicinity.  Another  point  that  should  not  be  overlooked  is  the 
probable  presence  of  an  appreciable  density  of  pairing  states  in  the 
electronic  structure  of  various  interfaces  such  as  oxide-semiconductor, 
metal-semiconductor  etc.  This  follows  since  these  interfacial  regions 
are  expected  on  the  whole  to  be  disordered  giving  rise  e.g.  to  weaker/ 
stronger  bonds.  Indeed  the  presence  of  such  centers  can  be  used  to 
understand  some  of  the  more  puzzling  electronic  behavior  of  the 
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localized  inversion  layer  regime  of  MOSFETs  where  one  is  dealing  with 
an  oxide-semiconductor  interface  in  contact  with  a quasi-two  dimensional 
electron  gas.  Furthermore,  recently^  we  have  carried  out  an  analysis 
of  the  origin  and  role  of  such  states  at  metal-semiconductor  (Schottky) 
interfaces  and  the  resultant  picture  has  been  found  to  be  consistent 
with  the  so-called  Covalent-Ionic  trend  . ^ Thus,  although  interface  or 
contact  effects  are  usually  ignored  we  expect  that  such  systems  should 
also  exhibit  a dielectric  loss  obeying  the  "universal  law"  and  a 
systematic  study  of  the  details  could  provide  a powerful  probe  of  the 
interfacial  structure. 

(b)  Dielectrics  with  Atom-atom  or  Molecule-molecule  or  Ion-ion 
or  Dipole-dipole  Interactions 

New  concepts  and  ideas  on  low-energy  excitations  in  real  glasses 

and  spin  glasses  have  been  recently  introduced  by  Anderson,  et  al,^ 

18  19 

Phillips  and  by  Anderson.  They  propose  the  existence  of  a 
statistical  distribution  of  localized  tunneling  levels  and/or  modes. 

A tunneling  mode  in  a real  glass  is  realized  by  an  atom  (or  group  of 
atoms)  which  has  an  energy  E(x)  as  a function  of  its  generalized  posi- 
tion  coordinate  x which  exhibits  two  local  minima  of  energy  difference 

v/' 

Z!E  separated  by  a barrier.  Similarly  in  spin  glasses  spins  are  con- 
sidered as  classical  dynamical  quantities  with  a potential  energy  sur- 
face that  is  a function  of  the  simultaneously  specified  orientations  of 
all  the  spins  (i.e.  a N-dimensional  configuration  space);  local  minima 
in  the  energy  correspond  to  metastable  states  of  the  spin  glass  associ- 
ated with  different  spin  configurations.  A tunneling  mode  for  spin 
glasses  ’ is  defined  in  spin  configuration  space  as  two  local  minima 
separated  by  a quantum-mechanical  energy  barrier.  Tunneling  between  one 
local  minimum  and  another,  if  it  occurs,  involves  the  rearrangement  of 
several  spins.  The  linear  specific  heat  observed  in  real  glasses  (spin 
glasses)  comes  from  tunneling  modes  whose  energy  barriers  are  suffi- 
ciently great  so  that  resonant  tunneling  of  atoms  (spins)  between  local 
minima  does  not  occur,  but  sufficiently  small  such  that  tunneling 
between  the  two  levels  can  take  place  during  the  time  span  of  the 
specific  heat  measurement.  Tunneling  modes  that  contribute  to  the  low 
temperature  linear  specific  heat  have  a density  of  levels  N(AE)  per  unit 
ZE  which  is  non-zero,  smooth  and  continuous  for  $ kT.  Those  tunnel- 
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ing  modes  that  contribute  to  the  low  temperature  linear  specific  heat 

17  19 

compose  only  a small  subset  ’ of  the  total  density  of  alternate 

17  19 

states  or  modes  with  level  splitting  Ae.  It  has  been  pointed  out  ’ 

that  there  are  also  a large  number  of  modes  having  small  AE  which  have 

their  two  alternate  states  inaccessible  to  each  other  because  their 

energy  barriers  are  too  large  for  tunneling  to  occur.  Those  pairs  of 

levels  are  practically  not  connected,  and  some  of  them  contribute  to  the 

zero  point  entropy  of  the  glass.  Indeed  experimental  measurements  of 
. . 20  21 

fused  silica  and  glycerol  has  shown  that  the  zero-point  entropy  is 
finite  for  both. 

Let  us  examine  the  transient  response  of  the  tunneling  modes  to 
sudden  potential  change  caused  by  fast  quantum  transition  of  some 
charged  species.  Tunneling  modes  whose  alternate  states  are  thermo- 
dynamically accessible  can  be  eliminated  at  the  outset  for  consideration 
of  infrared  divergent  response  at  oj/2tt  < 10  GHz.  Our  interest  is  in  the 
low  frequency  dielectric  response  where  co  is  smaller  or  much  smaller 
than  10  GHz  and  the  ambient  temperature  is  usually  room  temperature. 

Any  infrared  divergence  had  it  existed  would  be  obliterated  by  the 
effects  of  finite  temperature  T which  replaces  the  characteristic  t n 
by  an  exponentially  damped  dependence  in  the  response  function  ' for 
t »"R/kT.  Thus  we  need  to  consider  only  very  low  energy  tunneling 
modes  where  the  two  alternate  levels  of  a mode  must  be  thermodynamically 
inaccessible  to  each  other.  The  very  low  of  the  tunneling  modes 
guarantees  contribution  to  the  dielectric  response  at  corresponding  low 
frequencies  u)  s £E/1f,  and  thermodynamical  inaccessibility  enables  the 
infrared  divergence  to  survive  at  finite  T.  This  class  of  tunneling 
modes  should  exist.  Anderson  pointed  out  that  since  the  configura- 
tions of  the  atoms  (spins)  is  random,  there  must  be  very  many  locations 
(of  order  N,  the  number  of  atoms  or  sets  of  atoms)  where  there  are  two 
possible  configurations  of  very  similar  energies  E^  and  E2*  If  E^  and 
E2  are  independent  random,  variables,  then  the  probability  p(AE)  of 
finding  = ^-Ej  is  finite  as  £E  -►  0.  But  physically  this  is  not 
true  because  it  is  possible  to  tunnel  between  the  two  alternate  levels 
with  a tunneling  matrix  element  T^  even  though  it  is  small  for  thermo- 
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dynamically  inaccessible  tunneling  modes.  The  energy  level  separation 
will  be  at  least  AE>  | 2 1 » the  off-diagonal  matrix  element  between  the 
alternate  levels.  For  this  class  of  very  low  energy  inaccessible  tun- 
neling modes  (i.e.,  < 10  GHz)  the  physical  energy  difference  AE  is 
determined  by  the  off-diagonal  matrix  element  AE  = | 2 1 • 

Let  us  confine  our  further  discussions  to  only  spin  glasses.  There, 
it  has  been  argued  by  Anderson,  that  T^  being  a complex  matrix  element 
acts  like  the  x and  y components  of  the  random  field  that  prevents  the 
actual  level  splitting  Ae  going  to  zero  even  though  lE^-Ejl  ->-0  unless 
T12  * ® also*  Eor  low  frequency  dielectric  response,  we  are  particularly 
interested  in  the  AE  = | 2 1 * 0 limit.  T^  consists  of  two  random 
variables  since  it  has  real  and  imaginary  parts.  The  probability  that 
the  mode  energy  AE  lie  in  the  interval  | T|  and  | T|  + d|  T|  is  propor- 
tional to  | T| d| T| . Hence  the  density  of  states  of  inaccessible,  very 
low  energy,  tunneling  modes  J\|cae)  is  proportional  to  aE.  Now  the  sudden 
potential  change  that  induces  transitions  between  the  two  alternate 
levels  should  not  depend  on  AE.  Here  we  envisage  the  sudden  potential 
change  inducing  a virtual  transition  of  one  level  to  an  excited  state 
which  has  easy  access  to  the  other  level.  In  other  words  the  sudden 
potential  change  introduces  an  additional  transition  channel  that  allows 
thermodynamically  inaccessible  tunneling  modes  to  contribute  to  low 
frequency  dielectric  response.  Within  this  class  of  tunneling  modes, 
the  condition  for  infrared  divergence  | | 2|\J(AE)  = nAE  is  satisfied. 

The  universal  law  follows  for  spin-glasses.  A model  of  spin  glasses  in 

. 23  . 

the  Ising  model  formulation  has  recently  been  analyzed,  m which  the 

exchange  interactions  J. . are  assumed  to  be  distributed  randomly  and 

independently  of  one  another  over  both  positive  and  negative  values. 

Many  features  of  the  spin  glasses  including  the  low  temperature  linear 

specific  heat  have  been  derived  from  the  model. 

The  spin-glass  system  and  the  resultant  spin-spin  interaction 

models  can  often  be  transcribed  to  other  physical  models  with  non-spin 
. . 24 

interactions.  Well  known  examples  include  the  Ising  model  equivalence 
to  a lattice  gas  and  to  a binary  alloy.  A lattice  gas  is  a collection 
of  atoms  (molecules)  whose  positions  can  take  on  only  discrete  values 
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which  form  a lattice.  Each  lattice  site  can  be  occupied  by  at  most  one 
atom.  In  general  the  potential  energy  of  the  system  of  atoms  cor- 
responds to  a gas  in  which  the  atoms  are  located  only  on  lattice  sites 
and  interact  through  a two-body  potential  v(  - 5>j|)»  The  cor- 
respondence between  the  lattice  gas  and  the  Ising  model  is  seen  by 
identifying  occupied  sites  to  up  spin  and  empty  sites  to  down  spin  and 

the  nearest  neighbor  atom-atom  interaction  t.  to  -4  J..,  with  J. . the 

AA  1 j 1 j 

Ising  interaction  between  spins.  A binary  alloy  in  a lattice  model 
corresponds  to  sites  occupied  by  A or  B atoms  (molecules).  Let  e^, 
eAB’  Sb  rePresent  the  interaction  energies  between  the  atoms.  A site 
occupied  by  an  atom  A is  identified  with  an  up  spin  and  a site  occupied 
by  an  atom  B with  a down  spin.  The  quantity  (Ze^-e^-e^g)/ 4 then  cor- 
respond to  J in  the  Ising  model. 

Consider  dielectrics  where  atom-atom,  molecule-molecule  or  ion-ion 
interactions  are  important.  In  the  lattice  gas  and/or  binary  alloy 
modelling  of  dielectrics  with  random  interactions,  the  equivalence  to 
the  spin  glass  Ising  model  implies  a dielectric  state  corresponding  to 
the  spin  glass  state  exists.  Such  dielectrics  will  have,  in  analogy  to 
spin  glasses,  tunneling  modes  which  can  be  either  accessible  or  inac- 
cessible. In  direct  analogy  to  a tunneling  mode  in  spin  glasses  which 
corresponds  to  several  spins  turned  over,  in  these  dielectrics  a tunnel- 
ing mode  corresponds  to  the  change  of  the  atomic  (molecular  or  ionic) 
occupancy  of  several  sites  to  get  from  one  energy  minimum  to  the  other. 
The  essential  point  is  the  existence  of  very  low  energy  tunneling  modes 
in  these  dielectrics  which  are  only  accessible  when  a sudden  potential 
change  has  occurred.  This  class  of  tunneling  modes  again  satisfies  the 
criterion  for  infrared  divergence  and  hence  yields  the  universal  law. 

The  lattice  gas  and  binary  alloy  model  should  be  good  representations 

. . . . . . 25 

of  many  dielectrics  including  the  class  of  solid  state  ionic  conductor 

or  solid  electrolytes  such  as  Agl,  CaF  and  Na  6-alumina.  In  fact  ionic 

conductivity  for  these  solids  has  been  calculated  in  the  lattice  gas 

model.  In  the  case  of  Na  6-alumina,  there  is  the  repulsive  interaction 

among  the  diffusing  sodium  ions  and  also  the  attractive  interactions 

between  the  ions  and  their  randomly  distributed,  compensating  defects. 
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These  properties  imply  a lattice  gas  with  random  interactions.  There 

27  28 

is  indeed  ample  experimental  evidence  ’ for  the  existence  of  tunnel- 
ing modes  in  alkali  0-alumina  as  well  as  Ag  g-alumina.  In  particular 

28 

there  is  an  excess  low  temperature  specific  heat  contribution  which 
is  nearly  linear  in  T as  in  the  case  of  spin  glasses. 

To  conclude  this  Section,  we  note  that  the  apparent  arbitrary 
division  of  dielectrics  according  to  whether  electron  pairing  inter- 
actions or  ion-ion  interactions,  etc.,  dominate  the  behavior  of  the 
dielectrics  is  very  natural  after  all.  Ions  have  closed  atomic  shells 
and  molecules  are  usually  covalently  bonded.  In  both  cases  electron 
pairing  interaction  has  already  gone  to  completion,  although  the  origins 
of  the  pairing  interaction  in  the  two  cases  are  entirely  different.  The 
residual  interactions  are  then  the  ior-ion  or  the  molecule-molecule 
interactions,  which  then  should  play  the  important  role  in  providing 
correlated  states  and  their  excitations. 

IV.  Low  Frequency  Infrared  Divergent  Dielectric  Response 

Having  demonstrated  that  dielectrics  with  diverse  interaction  types 
should  have  invariably  some  very  low  frequency  excitations  that  con- 
tribute a time  dependence  of  the  form  t n at  large  t to  some  correlation 

function,  we  embark  on  the  derivation  of  the  dielectric  response  func- 

29  .... 

tion  and  examine  its  properties.  The  total  dielectric  polarization 

£ induced  by  an  electric  field  I£,(t)  can  be  calculated  by  standard 
methods^’^  of  linear  response.  The  interaction  of  the  polarization 
with  the  electric  field  is  given  by  = -£  . gL(t)  where  P is  the 

operator  of  the  polarizatiog.  The  perturbation  induces  a polariza- 
tion density  <P>=  <P>  + f 'i'(t-t')  . E (t')dt'  where  \b(t-t')  = 

-«P(t)  P(t')»  is  the  dielectfric  polarizability  tensor,  and  <P>„  is 
the  polarization  density  in  the  equilibrium  state  as  E > 0,  which  can 
be  nonzero  for  some  dielectrics  such  as  ferroelectrics . For  simplicity 
consider  the  dielectric  tensor  ^ to  be  diagonal.  In  the  case  when 
classical  statistical  mechanics^suf f ice  (as  often  is  the  case  for  die- 
lectrics at  finite  temperatures),  the  response  function  simplifies  to 
the  time  correlation  function  'f. . ( t— t * ) = 0 < P.(t)  where 
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<...>  denotes  averaging  with  the  equilibrium  distribution  function, 

® • 28 
0 ■ 1/k  T and  P.(t')  the  derivative  of  P.(t')  with  respect  to  t'. 

D 1 1 

If  P^(t)  takes  on  either  of  two  values  + pQ  and  makes  transitions 
from  one  value  to  the  other,  as  in  the  case  of  a system  of  particles  with 
a dipole  moment  or  the  case  of  a charged  particle  that  can  occupy  one 

of  two  alternate  sites,  then  ip. . can  be  readily  calculated  by  general- 

31  11  . . . . 

izing  the  method  to  take  into  account  a time  dependent  jump  transition 

rate  W(t).  Rewriting  t-t'  as  T,  we  wish  to  calculate  >p..(T)  = 

- (3<g^(t  g^(t-T)>o,  where  the  derivative  is  now  with  respect  to  X. 

Doing  this  we  obtain  the  result  'l^Cx)  =*  2f3p^  W(T)  exp(-2QJTW(X)dX)  for 
the  time  dependence  of  the  dielectric  repsonse  function.  The  task  that 
remains  is  to  calculate  W(x)  including  the  possibility  of  an  infrared 
divergence  of  correlated  states  excitations.  Let  <J*(x)  describe  the 
time  response  of  the  correlated  states  to  the  sudden  jump  of  the  elec- 
tron (dipole)  from  one  position  to  another  with  probability  per  unit 
time  Wq.  The  form  of  4>( t)  has  been  given  ’ and  in  our  notation  is 
<p(t)  * J®c  V^Jfc)  ( 1-e  ^C)dE/E^.  We  have  seen  in  the  last  section 
that  there  exists  some  class  of  correlated  states  in  the  dielectrics  we 
considered  so  that  V^Jfe)  =bV^E  is  proportional  to  E and  satisfies  the 
condition  for  infrared  divergence  in  the  number  of  these  low  energy  cor- 
related states  excitations.  The  integral,  <P(X),  has  been  evaluated 

under  this  circumstance  and  yields  <t>(x)  = bV  \Y+  ln(iE  t)  + 

J o c 

E,(iE  t)}  where  y = 0.5722,  E.(ix)  is  a standard  integral  which  vanishes 
at  large  x.  The  jump  transition  rate  is  W(x)  = WQ|e  |.  On 

defining  a time  Tq  by  1/X  = 2Wq  and  combining  equations,  we  obtain 

4<i i ( X ) 3 ( Bp^/ XQ)  |exp(— 4>(  x))  | exp(-/T|exp(-<i>(x))|dx/Xo)  (5) 

Consider  the  case  when  either  the  infrared  divergent  correlated  states 

2 

do  not  exist  or  the  coupling  VQ  of  the  hopping  charges  (dipoles)  to  the 

correlated  states  is  vanishingly  small.  Then  in  either  case  ^(x)  0 

2 

and  0.  • ( x)  * 3p  /x  exp(-  x/x_;  whose  Fourier  transform  is  Y..(u))  « 

2 'll  o o o t.1 

gp^(l+iojX0)  which  is  the  classical  Debye  susceptibility.  Rec  pturing 
the  classical  Debye  laws  by  turning  off  the  low  energy  correlated  state 
excitation  is  of  course  no  surprise.  The  interesting  point  is  that 
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dielectrics  or  dielectric  interfaces  in  nature  seldom  obey  the  Debye  law 
which  implies  there  should  exist  some  low  energy  correlated  states  exci- 
tations which  are  coupled  to  the  carriers/charges/dipoles  of  the  die- 
lectric. 

The  dielectric  response  function  for  E£t  »1  is  vJk^t)  = 

(6p2/t  )e  (E  t)  n exp(-e  n/(l-n)  T E n)  where  we  have  put  n = 

oo  o c 

bVQ  and  assumed  n < 1.  By  inspection  one  can  observe  that  although  the 
(Ect)  n term  may  initially  determine  the  T-dependence  of  \Jk^,  for  suf- 
ficiently large  values  of  T be  dominated  by  the  exponential  func- 
tion. This  occurs  roughly  at  T ^f(l-n)  en^  En  T n.  ¥•  •(<*))>  the 

PL  ^ ® J f i 

Fourier  transform  of  ^^(t)  of  Eq.  (5)  can  be  obtained  numerically. 
Several  representative  results  for  reasonable  choices  of  the  parameters 
n,  E and  T are  shown  in  Fig.  5.  A peak  in  x"(lo)  exists  and  its  loca- 


tion is  close  to  the  value  of  to  = 1/t 

P P 


1-n 


This  post-l/u)  “ peak  may  be 
identified  with  the  a or  the  8 peaks  commonly  observed  in  dipole  systems 
such  as  polymers,  liquids,  p-n  junctions,  ferroelectrics , liquid 
crystals,  cryogenic  polymers  and  some  glasses.  The  approximate  peak 
position  cSp  = £(l-n)en^  ToE^j^^n  ^ i-3  a decreasing  function  of 
increasing  tq  and  Ec  and  depends  sensitively  also  on  the  infrared 
divergence  exponent  n.  In  general  Tq  is  temperature  dependent  and 
usually  has  a clearly  defined  activation  energy  EA:  Tq(t)  = 


exp(EA/kgT) . This  alone  introduces  a temperature  dependence  into 

*/*  • . * - . . 'S' 


U)p  « exp(-EA/(l-n)kgT) , with  an  apparent  activation  energy  EA  of 
EA/(l-n).  Increase  in  temperature  will  cause  a lateral  shift  of  the 
universal  law  and  its  post-peak  along  the  frequency  axis.  If  a post- 
peak either  occurs  at  too  low  a frequency  to  be  measured  or  obscured  by 
another  mechanism,  it  can  be  revealed  at  higher  frequencies  by  raising 
the  temperature.  This  behavior  has  been  seen  for  example  in  the  ionic 
conductor  Hollandite  of  the  composition  gMgQ  gTi^  as  T ranges 

from  77K  to  230K.  If  n is  close  to  unity,  1?A  can  be  much  larger  than 

E.  and  corresponds  to  an  unreasonably  high  "activation  energy"  as  is 

A 1 
found  for  the  a loss  peaks  of  many  materials. 

A wide  range  of  dielectrics  have  associated  with  them  the  presence 
of  charge  carriers  of  electronic  or  ionic  nature.  These  charge  carriers 
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are  also  evidently  responsible  for  dc  conductivity.  Thus  one  expects 
that  charge  carrier  hopping  transitions,  under  excitation  by  a time- 
varying  electric  field,  do  not  necessarily  involve  only  two  preferred 
sites.  Consider  the  charge  carriers  that  do  not  jump  randomly  between 
two  states/sites,  then  the  dielectric  loss  is  simply  proportional  to 
the  probability  of  exciting  low  energy  correlated  state  excitations. 

With  the  same  time  response  function  of  the  correlated  states  <j>(t)  as 
displayed  in  preceding  paragraphs,  x"^)®  dt  exp(iTU))  exp(-<t>(l)) . 

For  EcT  large,  <J>(x)  can  be  approximated  by  ny  + n®-n(iE  X).  The  approxi- 
mate dielectric  loss  x"^)  is  then  proportional  to  l/u)^  n which  is 
identical  to  the  universal  law^  and  the  absence  of  a loss  peak.  This 
predicted  type  of  dielectric  response  is  indeed  observed  in  a very  wide 
range  of  dielectrics  of  all  physical  and  chemical  characteristics,  and 
interestingly  they  are  always  associated  with  the  presence  of  hopping 

charge  carriers  (Fig.  1).  A second  universal  law  (uj/u)  )n2  ^ will  follow 

n X ^ 

a first  (oj/co  ) 1 on  decreasing  u)  if  there  are  available  two  types  of 

correlated  states  that  can  contribute  to  infrared  divergences.  From  sum 

rule  considerations  on  x"^)*  we  expect  n2<n^  which  is  also  observed 

(Fig.  1). 

V.  SUMMARY  AND  DISCUSSIONS 

In  this  work  we  have  broadly  and  arbitrarily  classified  dielectrics 
according  to  the  type  of  interaction  or  correlations  inherent  in  all 
materials.  We  have  found  that  independent  of  the  type  of  correlations, 
a dielectric  in  general  has  gapless  "correlated  states"  whose  density 
of  states  is  continuous.  These  "correlated  states"  have  response  times 
much  longer  than  the  time  taken  by  the  hopping  between  sites  of  charged 
particles  or  jumping  between  orientations  of  dipoles.  Hence  the  hopping 
or  jumping  movements  can  be  considered  instantaneous  as  far  as  the 
"correlated  states"  are  concerned  and  they  experience  a sudden  change  of 
the  potential  induced  by  the  charged  particles  or  dipoles.  The  tran- 
sient response  of  the  system  is  the  emission  of  low-energy  excitations 
of  the  "correlated  states"  which  cause  the  response  to  have  a t n 
time  dependence  or  an  infrared  divergent  like  l/oo^  n frequency  response 
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of  Che  dielectric  loss.  We  have  thus  arrived  at  a fundamental  mechanism 
for  the  empirical  u)”  ^ dependence  (accompanied  sometimes  by  a peak  at 
low  enough  uj)  of  the  dielectric  loss  obeyed  by  nearly  all  dielectrics 
and  the  mechanism  is  operative  independent  of  the  type  of  physical 
structure  and  chemical  bonding  in  the  materials,  and  whether  the 
polarization  is  associated  with  permanent  dipoles  or  hopping  charge  car- 
riers of  electronic  or  ionic  nature. 

This  arbitrary  classification  of  dielectrics  according  to  the 
present  scheme  is  quite  general  after  all.  The  classification  is  based 
on  the  type  of  dominant  correlations  and  the  correlated  states  they 
render.  Detailed  developments  of  the  electron  pairing  correlations  and 
of  the  ion-ion  correlations  have  been  given.  Correlated  states  are 
identified  in  both  cases.  Types  of  correlations  other  than  those 
between  electrons  or  between  ions  could  conceivably  lead  to  some  sort  of 
"correlated  states"  as  has  been  demonstrated  explicitly  for  the  cases  of 
electron  pairing  correlations  and  the  ion-ion  correlation.  These  cor- 
related states  although  they  may  have  very  different  physical  origin  and 
interpretations  dependent  on  which  class  of  dielectrics  share  some  com- 
mon important  properties.  These  include  the  thermodynamical  inacces- 
sibility of  these  states  from  one  to  another,  and  yet  the  possibility  of 
excitation  when  a sudden  change  in  potential  occurs.  The  very  low 
energy  excitations  of  these  correlated  states  have  an  infrared  divergent 
behavior,  and  lead  to  the  low  frequency  dielectric  response  obeying  a 
universal  law,  x"(<i))  * 1/co^  n>  with  sometimes  the  appearance  of  a post 
peak  at  low  enough  u j.  The  Debye  law  holds  only  in  the  probably  seldom 

realized  cases  where  the  correlated  state  excitations  are  either  non- 
existent or  ineffective  because  of  weak  coupling  to  the  hopping  charges/ 
dipoles  that  contribute  to  the  Debye  susceptibility.  The  invariable 
deviation  from  the  Debye  laws  in  most  dielectrics  implies  that  the 
existence  of  very  low  energy  correlated  state  excitations  are  often  the 
rule  rather  than  the  exception.  We  emphasize  the  importance  here  of  not 
only  the  recognition  of  the  Curie-von  Schweidler  law  as  an  infrared 
divergence  phenomenon  but  also  the  subtle  task  of  identifying  the  (cor- 
related state)  excitations  that  are  responsible  for  it.  There  is  an 


21 


important  difference  between  the  present  case  and  the  Cerenkov  (or 
Bremstrahlung)  radiation  or  the  X-ray  edge  singularity  problem,  since 
energies  in  the  present  regime  of  interest  are  so  low  that  for  these 
cases,  the  spontaneous  photons  or  electron-hole  pairs  produced  infrared 
divergence  is  entirely  smeared  out  at  finite  temperatures.  This  is  not 
the  case  here  for  the  particular  correlated  states  responsible  for  such 
low  energy  dielectric  response  singularities  are  thermodynamically  inac- 
cessible from  one  to  another.  The  infrared  divergence  is  retained  at 
finite  temperatures  even  10  GHz.  In  all  infrared  divergence  problems, 
an  upper  cut  off  Ec  of  the  excitation  energies  E is  needed  either  to 
insure  convergence  at  large  E or  simply  that  we  run  out  of  these  excita- 
tions as  E increases.  In  the  specific  examples  we  discussed  here,  the 
upper  cut-off  energies  are  deduced  by  examining  the  nature  of  these 
excitations.  The  universal  law  x"(u))  “l/u)^  n may  be  modified  at  low 
enough  frequencies  in  dipolar  dielectrics  by  the  introduction  of  a peak, 
and  this  may  or  may  not  occur  within  the  frequency  spectrum  scanned, 
dependent  on  the  magnitude  of  Ec,  the  upper  cut-off  of  the  correlated 
state  excitations,  and  the  value  of  n.  The  occurrence  of  a post-peak  in 
some  classes  of  dielectrics  and  the  non-occurrence  in  other  classes  can 
be  correlated.  Order  of  magnitude  estimates  of  Ec  are  possible  for 
certain  classes  of  dielectrics  and  the  post-peak  frequency  predicted 
seems  to  be  consistent  with  experimental  data.  The  temperature  depen- 
dence of  the  post-peak  position  is  also  consistent  with  experimental 
data. 

In  addition  to  bulk  dielectrics  we  have  considered  also  the  inter- 
faces of  a dielectric  with  another  dielectric  or  a semiconductor  or  a 

metal.  An  interesting  example  is  the  thermal  oxidized  Si-Si0~  inter- 

. 13  . . 

face  in  MOS  device  structures.  The  present  authors  have  investigated 

the  local  electron  pairing  interaction  on  dangling  bonds  and  weaker/ 

stronger  bonds  (a  concept  also  introduced  by  Anderson^)  and  the 

. . 32 

resultant  electronic  structure  of  the  Si-SiC^  interface.  Both  the 
dangling  bonds  and  the  weaker/ stronger  bonds  give  rise  to  pair  states 
which  are  strongly  self-trapped  and  have  the  interesting  dynamic  charac- 
ter when  electrons  are  excited  in  pairs.  Correlated  pair  states  at  the 


22 


I 


interface  give  rise  to  electron  pair  excitations  with  arbitrary  low 
energies  and  should  give  rise  to  an  infrared  divergent  dielectric 
response.  We  wish  to  point  out  that  low  frequency  dielectric  response 
measurements  of  the  interfacial  region  could  be  a powerful  and  novel 
tool  for  the  characterization  of  devices.  These  measurements  may  have 
the  potential  of  yielding  more  in  depth  understanding  of  interfaces  when 
coupled  with  conventional  measurements  such  as  capacitance  versus  gate 
voltage . 
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Fig.  1.  A schematic  representation  of  the  various  observed  types  of 

dielectric  response  in  the  entire  range  of  solids.  The  upper 
set  of  diagrams  represent  the  shapes  of  the  logarithmic  plots 
of  x'  (co)-chain-dotted  lines,  and  Xn(a))”Solid  lines,  ranging 
from  the  ideal  Debye  through  the  Ct  and  8 peaks  and  on  to  the 
universal  dependence  for  charged  carrier  systems.  The  limiting 
forms  of  behaviour  are  represented  by  the  strong  low-frequency 
dispersion  with  small  values  of  n and  by  the  limiting  case  of 

j' 

frequency-independent  "lattice  response"  with  n ^ 1 . The  lower 
set  of  diagrams  represent  the  corresponding  complex  X plots. 

The  various  types  of  materials  obeying  the  respective  types  of 
response  are  shown  and  the  presumed  polarization  mechanisms  are 
indicated . 
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Fig.  5.  The  behavior  of  X' > X'  in  the  present  theory  for  several  dif- 
ferent values  of  n.  Note  the  peak  shape  is  independent  of 
a = e n^/  (l-njT^E^,11  but  strongly  dependent  on  n.  The  slope 
m of  each  of  these  log  (X1 ) versus  log  ( u>)  plots  varies  con- 
tinuously from  zero  to  one  for  log  (w)  < log  (o£ ) , where 
is  the  post  peak  position,  m for  a fixed  decrement  of  log  (^), 
i.e.  at  a value  of  (i)  with  log  (ot/ajp)  < 0 and  fixed)  decreases  as 
n increases.  In  view  of  this  property  one  should  not  take  the 
asymptote  of  the  lowest  available  frequency  measurements  of 
y"  (a))  and  attach  a universal  meaning  to  the  slope  of  that 
asymptote  but  rather  analyze  the  local  slope  m at  a fixed 
decrement  below  the  post  peak. 
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